Bose-Einstein Condensation of 88 Sr Through Sympathetic Cooling with 87 Sr 
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We report Bose-Einstein condensation of 88 Sr, which has a small, negative s-wave scattering 
length (ass = — 2ao). We overcome the poor evaporative cooling characteristics of this isotope 
by sympathetic cooling with 87 Sr atoms. 87 Sr is effective in this role in spite of the fact that it 
is a fermion because of the large ground state degeneracy arising from a nuclear spin of / = 9/2, 
which reduces the impact of Pauli blocking of collisions. We observe a limited number of atoms in 



the condensate (N n 
parameters. 



10 ) that is consistent with the value of ass and the optical dipole trap 



Bose-Einstein condensation of 88 Sr has been pursued 
for over a decade because of the promise of efficient laser- 
cooling to high phase space density using the (5s 2 ) -^o- 
(5s5p) 3 P\ narrow intercombination line [l| and loading 
of optical dipole traps that operate at the magic wave- 
length for this transition [2]. Recent interest in 88 Sr 
has focused on long-coherence time interferometers Q, 
optical frequency standards 0, Hj , and the existence of 
low- loss optical Feshbach resonances p, [7j. There has 
also been great interest generally in quantum degener- 
ate gases of alkaline-earth metal atoms and atoms with 
similar electronic structure because of potential applica- 
tions in quantum computing in optical lattices [8|4lQ| and 
creation of novel quantum fluids [ll[ • 

Early attempts to evaporatively cool 88 Sr to quantum 
degeneracy in an optical dipole trap [2|, LL2| were not suc- 
cessful in spite of initial phase space densities as high 
as 10 _1 , presumably because of a small elastic scatter- 
ing cross section. This was confirmed by measurements 
of the scattering lengths of all strontium isoto pes using 
photoassociative fl3l.ll3| and Fourier- transform |l5| spec- 
troscopy of Sr2 molecular potentials, which found that 
ass = — 2ao, where ao = 5.29 x 10 -11 m is the Bohr ra- 
dius. Here, we report Bose-Einstein condensation (BEC) 
of 88 Sr through sympathetic cooling with 87 Sr. 

Divalent atoms such as strontium and ytterbium [lH, 
17j often possess a large number of stable isotopes, which 
enables mass tuning of the s-wave scattering length. For 
strontium, the stable isotopes and abundances are 88 Sr 
(82.6%), 87 Sr (7.0%), 86 Sr (9.9%), and 84 Sr (0.6%). In 
such systems, the likelihood of finding an isotope with a 
scattering length that enables efficient evaporative cool- 
ing is very high, as was recently demonstrated through 
the condensation of 84 Sr (a 84 = 123 a ) 0, EH. For an 
isotope that has a poor scattering length for evaporative 
cooling, there are also numerous opportunities to find an- 
other isotope that can be used effectively for sympathetic 
cooling. For 88 Sr, the fermionic isotope 87 Sr is well-suited 
for this role. It has a large and positive s-wave scatter- 
ing length of as7 = 96 ao 13, IS] which leads to efficient 
thermalization and evaporation as long as the system is 
not highly polarized. The inter-isotope scattering length 
is also reasonable, as8-87 = 55 ao 3, EH, so that m a 
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FIG. 1: (color online) Partial level diagram for 88 Sr (- -) and 
87 Sr ( — ) showing hyperfine structure and isotope shifts [22I - 
[241 ]. Total quantum number F is indicated for 87 Sr levels. 



mixture, 88 Sr will be efficiently cooled (and evaporated) 
through collisions with 87 Sr. 

Bosons are normally used to sympathetically cool 
fermions [io| rather than the other way around, because 
identical fermions suffer from collisional Pauli blocking, 
which reduces evaporation efficiency in the quantum de- 
generate regime [21|. We do not observe significant limi- 
tations due to Pauli blocking in the experiments reported 
here, and we suspect this is because 87 Sr has a large nu- 
clear spin (/ = 9/2) and ground state degeneracy. This 
suppresses the Fermi temperature and allows 88 Sr to be 
cooled to high phase space density before Pauli blocking 
of 87 Sr collisions becomes important. 

Details about our apparatus can be found in refs. 
[14|, ll9|, |2S|. Formation of ultracold mixtures of stron- 
tium isotopes benefits from the ability to magneti cally 
trap atoms in the metastable (bsbp) 3 P2 state [12l l26M28| . 
which has a 10 min lifetime [29|. One isotope is trapped 
from a Zeeman slowed beam in a magneto-optical trap 
(MOT) operating on the (5 s 2 ) 1 So-fisbp) 1 Pi transition at 
461 nm. This transition is not closed and approximately 
1 in 10 5 excitations results in an atom decaying through 
the (5s4:d) 1 D2 state to the (5s5p) 3 P2 state, where it can 
be trapped in the quadrupole magnetic field of the MOT. 
After accumulating a desired number of atoms, limited by 
the loading rate and observed Z P^ lifetime of about 25 s, 
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the cooling laser frequency is then switched to cool and 
accumulate another isotope. In our experiment, we load 
88 Sr for 3 s and then 87 Sr [22| for 30 s, which yields an ap- 
proximately equal number of atoms of each isotope dur- 
ing evaporative cooling. The laser parameters for traj 
ping 88 Sr are given in ref. [26]. For trapping 87 Sr 
the laser is approximately 70 MHz red-detuned from the 
1 S (F = 9/2)- 1 Pi(F = 11/2) transition (slightly more 
than 2r, where T = 30.5 MHz is the natural linewidth of 
the transition [30| ) . 

3 P2 atoms are returned to the ground state using 60 ms 
of 3 W/cm 2 of excitation on the (5s5p) 3 P 2 -(5s4:d) 3 D2 
transition at 3 /xm [25|. The isotope shift (/g7 — /ss = 
110 MHz) [25] is small compared to the - 500 MHz width 
of the repumping efficiency curve [31] for 88 Sr and the 
~ 3 GHz width of the hyperfine structure in 87 Sr [25|. We 
tune the 3 jam laser 1.6 GHz blue detuned from the 88 Sr 
resonance, which optimizes the 87 Sr repumping while re- 
ducing the 88 Sr number by 80%. This is a reasonable 
compromise given that the number of 87 Sr atoms is the 
limiting factor in the experiment. The 461 nm MOT is 
left on at the optimal 87 Sr detuning to maximize the 
number of captured 87 Sr atoms, but this is only ~ 5 nat- 
ural linewidths red-detuned of the 88 Sr 1 6'o- 1 Pi transition 
[H, so it also serves to aid recapture of this isotope. 
We typically recapture approximately 1.1 x 10 7 88 Sr and 
3 x 10 7 87 Sr at temperatures of a few millikelvin. 

The 461 nm light is then extinguished and 689 nm light 
is applied to drive the (5s 2 ) 1 So-(5s5p) 3 P\ transitions and 
create intercombination-line MOTs for each isotope. The 
resonance frequencies in each isotope are well-resolved 
compared to the 7.4 kHz transition linewidth, so the si- 
multaneous MOTs are compatible with each other. The 
parameters of the 1 So- 3 -Pl lasers for 88 Sr and 87 Sr [31] 
are similar to the conditions in ref. [1] and ref. 
respectively. As many as 70% of the atoms are ini- 
tially captured in the intercombination-line MOT. (For 
experiments with one isotope, the loading phase and 
intercombination-line lasers for the other isotope are 
omitted.) 

After 400 ms of 1 So~ 3 Pi laser cooling, an optical dipole 
trap (ODT) consisting of two crossed beams is overlapped 
for 100 ms with the intercombination-line MOT with 
modest power (3.9 W) per beam. The ODT is formed by 
a single beam derived from a 20 W multimode, 1.06 jam 
fiber laser that is recycled through the chamber to pro- 
duce a trap with equipotentials that are nearly oblate 
spheroids, with the tight axis close to vertical. Each 
beam has a waist of approximately 90 jam in the trapping 
region. 

Immediately after extinction of the 689 nm light, the 
ODT power is ramped in 30 ms to 7.5 W to obtain a 
trap depth of 25 joK. Typically the atom number, tem- 
perature, and peak density at this point for both 88 Sr 
and 87 Sr are 3 x 10 6 , 7/iK, and 2.5 x 10 13 cm- 3 . The 
peak phase space density (PSD) for 88 Sr is 10 -2 . 
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FIG. 2: (color online) (A) Temperature evolution in an ODT 
with trap depth of U/kB = 23 /oK for samples of 88 Sr and 87 Sr 
alone and for each in a mixture. The number of each isotope 
present initially is approximately 10 6 . (B) Number and (C) 
temperature for a mixture along a typical forced evaporation 
trajectory. 



For diagnostics, we record 1 So- 1 Pi resonant absorption 
images of samples after a time of flight varying from 10 
to 40 ms. Because of broadening of the resonance due to 
hyperfine structure, 87 Sr atoms present would contribute 
significantly to the absorption when imaging at the 88 Sr 
resonance frequency. To remove 87 Sr atoms and obtain 
clean 88 Sr images, light resonant with the 1 So(F — 9/2)- 
3 Pi(F — 11/2) transition in 87 Sr is applied during the 
first 2 ms of the time of flight. 87 Sr atoms are imaged with 
linearly polarized light resonant with the 1 So(F = 9/2)- 
1 Pi(F = 11/2) transition, and the contamination due to 
88 Sr is small and easily accounted for (3l| . 

To investigate the collisional properties of the different 
isotopes and the mixture, the evolution of number and 
temperature were recorded in a fixed potential (Fig. [2^). 
For 88 Sr alone, evaporation is inefficient and a typical ra- 
tio of the trap depth to the sample temperature is rj « 4, 
as observed previously [12|. 87 Sr, however, approaches 
rj « 9. Modeling [33] of the free-evaporation trajectory 
for 87 Sr alone suggests a moderate degree of polarization 
[3l| that will be investigated in future studies. The tem- 
peratures of 87 Sr and 88 Sr atoms in a mixture with peak 
densities of 8 x 10 12 cm -3 track each other closely and 
approach rj « 8, indicating that 87 Sr provides efficient 
sympathetic cooling of 88 Sr. 

Figure [2] shows the number (B) and temperature (C) 
for a typical forced evaporation trajectory with a mix- 
ture. We decrease the laser power according to P = 
P /(l + t/T)P + P o ff set , with time denoted by t, /3 = 1.4, 
and r = 1.5 s. This trajectory without P Q ff set was 
designed [34] to yield efficient evaporation when grav- 



3 



ity can be neglected. Gravity is a significant effect in 
this trap for Sr, and to avoid decreasing the potential 
depth too quickly at the end of the evaporation, we 
set Poffset = 0.7 W, which corresponds to the power at 
which gravity causes the trap depth to be close to zero. 
The lifetime of atoms in the ODT is 30 s. This allows 
efficient evaporation and an increase of PSD by a factor 
of 100 for a loss of one order of magnitude in the number 
of atoms. The 87 Sr and 88 Sr remain in equilibrium with 
each other during the evaporation. 87 Sr atoms are lost 
at a slightly faster rate, as expected because essentially 
every collision involves an 87 Sr atom. 

Figure [3] shows false color 2-dimensional renderings of 
(left) and 1-dimensional slices through (right) the time- 
of-flight absorption images recorded after 16 ms or 22 ms 
of expansion for various points along the evaporation tra- 
jectory. At 5 s of evaporation, the distribution is fit well 
by a Boltzmann distribution, but at 6 s, a Boltzmann 
distribution fit to the high velocity wings clearly under- 
estimates the number of atoms at low velocity. A fit using 
the Bose-Einstein distribution [35], however, matches the 
distribution well. The fugacity obtained from this fit is 
1.0, indicating this is close to the critical temperature for 
condensation. With further evaporation, the presence of 
a Bose-Einstein condensate is indicated by the emergence 
of a narrow peak at low velocity, and eventually, a pure 
condensate is observed. 

At the transition temperature, 2 x 10 5 87 Sr atoms re- 
main at a temperature of 0.2 /aK. This corresponds to 
T/Tp = 0.9 for an unpolarized sample, which is non- 
degenerate and above the point at which Pauli blocking 
significantly impedes evaporation efficiency [2l| . 

88 Sr has a negative scattering length, so one expects a 
collapse of the condensate when the system approaches 
a critical number of condensed atoms given by [36| 



N cr = 0.575 



a>ho 



(i) 



Here a^ = [%/ (mcJ)] 1 / 2 is the harmonic oscillator length, 
where m is the atom mass, Ti is the reduced Planck con- 
stant, and uj = (uJxUJyUJz) 1 / 3 is the geometric average of 
the oscillator frequencies. One should also see large fluc- 
tuations in the number of condensed atoms during the 
evaporation due to repeated collapses and refilling of the 
condensate. To investigate this, we recorded the conden- 
sate number for various points in the evaporation trajec- 
tory over many experimental runs. For absorption images 
with a condensate and thermal pedastal, we fit the wings 
of the thermal cloud, which are beyond the condensate 
radii of about 23 /im, to a Bose distribution with fugacity 
set to 1. The residuals of the fit represent the conden- 
sate atoms which are fit with the standard Thomas- Fermi 
functional form [19|, |37j to determine their number. 

Figure H] shows the observed condensate number along 
the evaporation trajectory from 7 to 10 seconds, as well 
as maximum values N cr predicted by Eq. [1] where uj is 
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FIG. 3: (color online) Appearance of Bose-Einstein conden- 
sation in absorption images (left) and areal density profiles 
(right). Data correspond to 16ms (bottom) or 22ms (top 
three) of free expansion after indicated evaporation times (t). 
Images on the left have the same time stamp as on the right. 
The areal density profiles are from a vertical cut through the 
center of the atom cloud, and temperatures are extracted from 
2D Bose-Einstein distribution fits to the thermal pedestal. At 
7.5 s, a bimodal distribution is evident and indicative of Bose- 
Einstein condensation. For bimodal data, the central region 
is excluded from the fits and the fugacity is constrained to 1. 
A pure condensate is shown at 9 s of evaporation. 



14000 

12000 

LU 

c 10000 

g 8000 
< 

•5 6000 

1 4000 
z 2000 











▲ 


* Data 






— Model 






Model Confidence Interval 




± ▲ 


▲ 














% 


▲ ▲ 

... A. ... : 


▲ 










▲ 



1.2 1.3 1.4 

ODT Power [W] 



1.5 



FIG. 4: (color online) Comparison of observed condensate 
number and maximum condensate number predicted by a 
model of the trapping potential along the evaporation tra- 
jectory. 
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determined from knowledge of the ODT potential, with 
confidence intervals reflecting uncertainties of the waists 
of ODT beams, of about 10%, and the uncertainties of 
the scattering length a$s = — 2.0(3) ao [15]. As expected 
for attractive interactions, large fluctuations in number 
are observed. While the best guess curve for N cr falls 
below some of the data, the upper bound is reasonably 
well accommodated by the confidence interval. Uncer- 
tainties in knowledge of the trap become larger at low 
ODT power because of the increasingly important role of 
gravity, which weakens the trap. 

We have described the Bose-Einstein condensation of 
88 Sr through sympathetic cooling with 87 Sr. Observation 
of large fluctuations in number below a maximum number 
of condensed atoms is consistent with the small, negative 
value of <288- Because of the very weak interactions, it 
should be possible to change the sign of the scattering 
length with an optical Feshbach resonance 0, 0] while 
keeping induced inelastic losses low. This suggests many 
possible future experiments, such as creation of matter- 
wave solitons in two dimensions [Hj and quantum fluids 
with random nonlinear interactions (39| . 
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